HUB-EP-97/25 
HD-THEP-97/14 



The Two-Loop Euler-Heisenberg Lagrangian 
in Dimensional Renormalization 



D. Fliegner Q, M.Reuter ^, M. G. Schmidt f\, C. Schubert 

^Institut fiir Theoretische Physik, Universitdt Heidelberg 
Philosophenweg 16, D-69120 Heidelberg, Germany 

^ Deutsches Elektronen- Synchrotron DESY, 
Notkestrasse 85, D-22603 Hamburg, Germany 

'^Institut fiir Physik, Humboldt Universitdt zu Berlin 
Invalidenstr. 110, D- 101 15 Berlin, Germany 



Abstract 

We clarify a discrepancy between two previous calculations of the two-loop QED 
Euler-Heisenberg Lagrangian, both performed in proper-time regularization, by 
calculating this quantity in dimensional regularization. 
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One of the earliest results in quantum electrodynamics was Euler and Heisenberg's |T|] 
calculation of the one-loop effective Lagrangian induced by an electron loop for an electro- 
magnetic background field with constant field strength tensor F^^. Written in Schwinger's 
proper-time representation 0, this Lagrangian reads 

^(1) r^T L '^^r-m^T {eaT){ebT) 

^spin[ J g^2_/Q ^3^ tan(eaT) tanh(e6r) ■ ^' 
Here a and h may be expressed in terms of the two invariants of the electromagnetic field, 
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- b2 + a/(E2 - B2)2 + 4(E • B; 



h 
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(E2 - b2) + ^(£2 - B2)2 + 4(E • B)^ 



(2) 



Schwinger also supplied the corresponding result for scalar quantum electrodynamics pi, 



m = J_ r {eaT){ebT) 
'^^^i^ ^ 167r2io T3 sin(ear) sinh(e6r) ■ 

Even the first radiative correction to the Euler-Heisenberg Lagrangian, due to the exchange 
of one internal photon in the electron loop, has been obtained many years ago by Ritus |^ 
in terms of a two-parameter integral. 

Dittrich and one of the authors |Q later obtained a similar but simpler integral represen- 
tation for the same quantity. They also verified the agreement of both representations in the 
strong field limit. 

In a recent publication |5| three of the present authors showed that this type of calculation 
can be considerably simplified using the worldline path integral variant j^, |^, |8|, |ll], 12 1 



of the Bern-Kosower formalism [13, 14]. 



This third calculation led to exactly the same parameter integral for the regularized ef- 
fective Lagrangian as the Dittrich-Reuter calculation, however in a more elegant way. 

All three calculations were performed in four dimensions, proper-time regularization, and 
on-shell renormalization. This choice of regularization keeps the integrations simple, but at 
the two-loop level makes it already somewhat nontrivial to achieve a consistent on-shell renor- 
malization. The difficulty is due to the non-universal nature of the proper-time cutoff, and 
somewhat similar to the problems encountered in multiloop Feynman diagrams calculations 
performed with a naive momentum space cutoff (see [^] and references therein). 

In fact, the calculation in |^ was incomplete in so far as we were not able there to deter- 
mine the finite part of the one-loop mass displacement appropriate to the present calculational 
scheme. Moreover, as a byproduct of this investigation we found that the formulas obtained 
by Ritus and Dittrich-Reuter for the renormalized effective Lagrangians are incompatible as 
they stand, and for the very same reason; if at all, they can be identified only after a certain 
finite mass renormalization. One had to conclude that in (at least) one of the two previous 
calculations the physical renormalized electron mass had been misidentified (the strong field 
limit checked in Q is not sensitive to this discrepancy). 

In the present letter, we clarify this matter by recalculating this effective Lagrangian in 
dimensional regularization. This paper should thus be seen as supplementing |^, and the 
reader is referred to that publication for some of the details of the formalism used here. 
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As usual in applications of the "string-inspired" technique to QED ||T^, our calculation 
of the two-loop Euler-Heisenberg Lagrangian for spinor QED will yield the corresponding 
Lagrangian for scalar QED as a partial result. We therefore consider the scalar QED case 
first. 

In scalar QED, our starting point is the following worldline path integral representation 
for the bare dimensionally regularized two-loop Euler-Heisenberg Lagrangian [p|. 
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Here T denotes the usual Schwinger proper-time parameter for the scalar circulating in the 
loop. The path integral T>y is to be performed over the space of all closed loops in D - 
dimensional spacetime, with an arbitrary but fixed center of mass xq, and traversed in the 
fixed proper-time T. The parameters Ta^b parametrize the end-points of the photon exchanged 
in the scalar loop. This form of the photon insertion corresponds to Feynman gauge. We use 
euclidean conventions both on the worldline and in spacetime. 

This representation of the effective action in terms of a first-quantized path integral goes 
essentially back to Feynman ||l^, except that we have dimensionally continued it to D 
dimensions. 

It is useful to take the background field A in Fock-Schwinger gauge centered at the loop 
center of mass xq 0], where one has A^{xq + y) = \y^Fi,^. Exponentiating the denominator 
of the photon insertion using a proper-time parameter T, one arrives at 



ai^-m^T / dfiinfy- / dTa / dn 
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The new path integral is gaussian, so that its evaluation amounts to a single Wick contraction 
of {ija -yb)- This leads to ||] 



xdet 2 



2 Jo ^ 
s'm{eFT) 



2'r„ D 

T-- 



eFT 



det 2 



:C. 



OO _ rT rT 

dT / dTa / dn 
Jo Jo 



ab 



iija ■ i/b) 



(6) 



Of the (Lorentz) determinant factors appearing here, the first one yields just the one-loop 
Euler-Heisenberg-Schwinger integrand eq. (P). It represents the dependence of the path inte- 
gral determinant on the external field while the second determinant factor represents its 
dependence on the photon insertion. The Wick contraction is given by |@] 



iija ■ yb) = tr 



Qsab + 



1 {GBaa — Q Bab) {Q Bab ~ Ssbb) 

2 f — \Cab 



(7) 
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where 
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eF 



-ieFTGn,, + ieFGB 



(8) 



is the bosonic one- loop worldhne Green's function modified by the constant field fs, 17, 19 1, 
and Cab = Qsaa — Qsab — Qsba + QBbb- Qb generalizes the ordinary worldline Green's function 
Gb, 



GB{Ti,Tj 



Ti - Ti 



(9) 



We will often abbreviate Gbi2 = Gb{ti-,T2) etc., and a "dot" always denotes a derivative 
with respect to the first variable, e.g. Gbi2 = sign(ri — T2) — 2^^3.^^. 

Performing a partial integration with respect to on the first term in eq. (0) one can 
derive the alternative parameter integral 
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For our present purpose, we can restrict ourselves to the pure magnetic field case. This field 
we take along the z-axis, so that F^^ = B, F^^ = —B are the only non- vanishing components 
of the dimensionally continued field strength tensor. 
We also introduce the following abbreviations, 



z 
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eBT 

{f + Gsab) 
{f + G%ab) 



and the z - dependent Green's function G 
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(11) 
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(compare eq. (HD). 

With these definitions, we can then rewrite the various traces and determinants appearing 
in eqs. (P), (lTo|) as 
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(13) 



The term involving (5(T(j — r;,), stemming from t/Bafei can be omitted, since it will not contribute 
in dimensional regularization (it corresponds to a massless tadpole insertion in field theory). 

Inserting these expressions into either eq. (^ or eq. (|lO|), one finds that the resulting 
integrals suffer from two kinds of divergences: 

1 . An overall divergence of the scalar proper-time integral dT at the lower integration 
limit. 

2. Divergences of Jq dTa Jq dri^ at the point = where the photon end points become 
coincident. 

The first one must be removed by one- and two-loop photon wave function renormalization, 
the second one by the one-loop renormalization of the scalar mass. 

It turns out that this program is easier to carry out on a certain linear combination of 
eqs. (^) and (|lO|), namely 



^^jy^ X eq.(6) + j^-x eq.(0) 



(14) 



(A similar simplification can be achieved by taking the photon insertion in Landau gauge, 
though the resulting parameter integrals are not identical). Moreover, we rescale to the unit 
circle, Ta^b = Tua,b,T = TT, and use translation invariance in r to set = 0. Thus in the 
following we have G^afe = Ma(l - Wa), G^afe = 1 - 2Ma- 
The resulting integral we write 



dT 



-m^Trp2-D 



oo rl 

df / dUaI{z,Ua,f,D) (15) 
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where the rescaled integrand I{z, Ua,T, D) depends on T only through z. 

In contrast to the calculation in proper-time regularization, the T - integration is nontriv- 
ial in dimensional regularization. It will therefore be easier to extract all sub divergences before 
performing this integral. An analysis of the divergence structure shows that the integrand 
can be rewritten in the following way. 



K{z,Ua,D) 



I 1 1J_ cy IJ_ 

dfl{z,Ua,f,D)=Ko2{z,Ua,D) + f{z,D)G^J +0{z\G^J) (16) 
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Kq2 consists of the terms constant and quadratic in z, which are the only ones causing a 
divergence at T = 0. The second term is 0{z'^), so that its integral already converges at 
T = 0, however it diverges at = 0, 1. 

After splitting off these two terms, the integral over the remainder is ah'eady finite, so 
that one can set D = 4 in its exphcit computation. For D = 4 the T - integral becomes 
elementary, and yields 
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The divergences will now be removed by mass and photon wave function renormalization, 
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So far we have worked in the bare regularized theory, so that all our previous formulas should, 
for the following, be considered written in terms of mo,eo,-Bo instead of m,e,B (note that 
this leaves z unaffected). 

Since we aim at a direct comparison with previous calculations, the renormalization will 
be done using on-shell rather than minimal subtraction. In on-shell subtraction, the photon 
wave function renormalization has the effect of simply removing the z'^ - part of Kq2, and the 
remaining z - independent term can, of course, be also discarded. 

The removal of the divergence caused by the second term in eq. (|l6|) takes more effort, 
and the mere possibility requires a little conspiration. 

Let us denote the corresponding contribution to the effective Lagrangian by Gsca.i{z, D). 
For this term the Ua - integration factors out, yielding 
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G^J = 5 (2 - - , 2 - -) = -- + + 0(6) 



(20) 



/o — V 2' 2/ 

where B denotes the Euler Beta- function. 

To proceed further, it is essential to note that the function f{z,D) can be related to 
the integrand of the scalar one-loop Euler-Heisenberg Lagrangian, eq. (|3|). In dimensional 
regularization, and with the two terms lowest order in z subtracted out via one- loop photon 
wave function renormalization, this Lagrangian reads 
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On the other hand, we can rewrite 
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By a partial integration over T, we can therefore reexpress 
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(there are no boundary terms since f{z) = 0{z'^)). 

At the two-loop level, the effect of mass renormalization consists in the following shift 
produced by the one- loop mass displacement Sttiq, 



(24) 



Sttiq is generated by the UV divergence of the one- loop scalar self energy in scalar QED. This 
quantity we have to take from standard field theory. In dimensional regularization one has |^ 



X 2 2^0 



-- + 7 - 3[7 - ln(47r)] - 31n(m^ 
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Here e = D — 4, and 7 denotes the Euler-Mascheroni constant 0. Expanding eqs. (|2^), ( p3| 
and (|2^) in e one finds that, up to terms of order 0(e), 
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Here this denotes the mass displacement itself, 



^Note that this differs by a sign from Sm^ as used in 
there the corresponding counterterm. 

■^In comparing with M, U, fel bol note that there this constant had been denoted by ln(7). 



(26) 
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Note that the whole divergence of Gscai{z, D) for D 4 has now been absorbed into 6mQ. 

Our final answer for the two-loop contribution to the finite renormalized scalar QED 
Euler-Heisenberg thus becomes 
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(27) 



As far as is known to the present authors, the only previous calculation of the two-loop Euler- 



Heisenberg for scalar QED is the one in |20, 21|. The parameter integral given there is rather 
different from ours, and we have not succeeded in directly identifying both representations. 
However, we have used MAPLE to expand both formulas in a Taylor expansion in B up to 
order 0{B^^), and found exact agreement for the coefficients. 
Let us just give the first few terms in this expansion. 
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(28) 

The expansion parameter has been rewritten in terms of Bcr = ~ 4.4 • lO^^G. 

The corresponding calculation for spinor QED is completely analogous. In the worldline 
superfield formalism of |^, P, 22|, the formulas (^), (0) generalize to the following integral 
representation for the two- loop effective action induced by the spinor loop. 
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The Q appearing here is the constant field worldline superpropagator, 

O{ri,9i;T2,02) = Qb{ti,T2) + Oi020F{n,T2) 

which besides the bosonic propagator eq. (^) also contains a fermionic piece, 

Qf{ti,T2) = sign(ri - T2) 



(30) 



(31) 



cos (eFT) 

Our superfield conventions are D = — 9-^, J d99 = 1. 

Performing the Grassmann integrations, and removing Gsab by partial integration over 
Ta, we obtain the equivalent of eq. (|lO|) , 
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Note that this formula reduces to eq. (^), if one replaces tan(eFr) by sin(eFr), and deletes 
all the ^i?, as well as the global factor of —2 (which accounts for the difference in statistics 
and degrees of freedom between the spin and spin ^ loops). 

Contrary to the scalar QED case, here the partially integrated integral is already a suitable 
starting point for renormalization (for more on this point see chapter 7 of Q). 

Specializing to the magnetic field case, it is again easy to calculate the Lorentz determi- 
nants and traces. After rescaling to the unit circle, one obtains a parameter integral 
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The extraction of the subdivergences yields 
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Lo2 is again removed by photon wave function renormalization. Denoting the contribution 
of the second term by Gspin{z, D), we note that the Ua - integral is the same as in the scalar 
QED case, eq. (|20|). Using the following identity analogous to eq. (p^), 
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(36) 



we partially integrate the remaining integral over T. The ^ - part of Gspin is then again found 
to be just right for absorbing the shift induced by the one- loop mass displacement. 
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Up to terms of order e one obtains 
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Our final result for the on-shell renormalized two-loop spinor QED Euler-Heisenberg La- 
grangian is 
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Comparing with the previous results by Ritus and Dittrich-Reuter, we have again not suc- 
ceeded in a direct identification with the more complicated parameter integral given by Ritus 
j^]. However, as in the scalar QED case we have verified agreement between both formulas 
up to the order of 0{B^^) in the weak-field expansion in B. The first few coefficients are 
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On the other hand, our formula almost allows for a term by term identification with the 
result of Dittrich-Reuter as given in eqs. (7. 21), (7. 22) there. This requires a rotation to 
Minkowskian proper-time, T ^ is, a transformation of variables from Ua to v := GBabi the 
use of trigonometric identities, and another partial integration over T for the last two terms in 
eq. (^). The only discrepancy arises in the constant 18, which reads 10 in the Dittrich-Reuter 
formula. 

Since this constant can be adjusted by a change of the finite constant appearing in 6mo, we 
conclude that the two previous results for this effective Lagrangian differ precisely by a finite 
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mass renormalization. Moreover, it is clear that Ritus' formula is the one which correctly 
identifies the physical electron mass. 

Finally, let us mention that the one-loop Euler-Heisenberg Lagrangian, and perhaps even 
its two-loop correction considered here, may possibly be measured in optical experiments in 
the near future p^ ]. 
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